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, Abstract. We report an extension of the source imaging method for imag- 

ing full three-dimensional sources from three-dimensional like-pair correlations. 
Our technique consists of expanding the correlation data and the underlying 
' source function in spherical harmonics and inverting the resulting system of 

\Q , one-dimensional integral equations. With this method of attack, we can image 

the source function quickly, even with the extremely large data sets common 
in three-dimensional analyses. We apply our method to the recently measured 
E859 un-Coulomb corrected data. 
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1. Introduction 

A single nucleus-nucleus collision at RHIC can produce tens of thousands of parti- 
cles. Basic quantum mechanics tells us that all of these particles are entangled; they 
are all part of the same final state wavefunction! The entanglement is a result of 
conservation laws, particle statistics and inter-particle interactions. Intensity inter- 
ferometry uses pairs entangled by particle statistics and inter-particle interactions 
to measure the space-time extent of a nucleus-nucleus collision, i.e. the source. This 
technique starts with the construction of the correlation function: 

dNp a i TS I eOVmlxed 
C W = A^ A- / A^ A^ ■ 1 

dp\dp 2 I dpxdpi 

Here, the numerator is the spectrum of particle pairs from the same event and 
constitutes the signal as these pairs are entangled. The denominator is the spec- 
trum of particle pairs from different events, using the same cuts as the signal, and 
constitutes the background as it is constructed to look like the signal without the 
entangled pairs. 
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For identical charged mesons, the correlations are nearly Gaussian except near 
the origin where the Coulomb force distorts the correlation. This effect is usually 
weak and corrected out leaving a correlation one fits to a Gaussian: 

C(0 = l + Aexp(-J r - where ?=±(ft -ft). (2) 

Here the correlation is tabulated in the Bertsch-Pratt coordinates (namely the side- 
ward, outward and longitudinal directions). Following uncertainty principal-type 
arguments, the radii of the correlation are taken as the radii of the source. Unfor- 
tunately, the Coulomb correction has a tendency to distort the extracted radii [Q. 
Usually the distortion is unimportant, but at RHIC we find that the sideward and 
outward radii are so close that this distortion can make the difference between radii 
that seem to indicate imaginary freeze-out durations and radii that do not. 

For identical charged mesons, we would like to forget (or at least correct away) 
final-state interactions. For other pairs such as proton pairs or proton-A pairs, final- 
state interactions dominate the signal [ E] • The Koonin-Pratt equation provides 
more general framework [Ej: 

C% = J d 3 r\^(f)\ 2 S(r). (3) 

Here S(r) is the source function and it gives the probability of producing a pair a 
distance f apart in the pair CM frame. Also, here 3>i \r) is the pair's final state 
wavefunction which we can compute if we know their interaction. In this paper, we 
will refer to the square wavefunction as the kernel, K{q, r) = |$l '(r)| 2 . 

Since we can no longer read off the radii to get the source size, we must invert 
Eq. (J3J to get S(r) directly. We will do this with an extension of the one-dimensional 
source imaging technique introduced in Refs. [0]. Source imaging allows us to 
cleanly separate the effects due to final-state interactions and symmetrization from 
effects due to the source function itself. Because we make no assumptions about 
the functional form of the source function, the imaging can reveal non-Gaussian 
sources. In the remainder of this paper, we will outline the three-dimensional imag- 
ing process. We first state our test problem, then detail the imaging process and 
finally show the imaging results on the test problem. With this experience in hand, 
we analyze data from the AGS-E859 experiment. 

2. Setting Up the Problem 

We begin with an emission function that has finite lifetime (r — 20 fm/c ), finite 
temperature (T = 50 MeV), and a Gaussian spatial profile (R x , Ry,Rz = 4 fm): 

^'.a««pH/T)cx p (-5^) «p (-^-^-^). (4) 

1 We choose r to be in the ballpark of the ui meson lifetime (23 fm/c). 
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On the face of it, this source would seem to produce an unchallenging spherically 
symmetric source function. However since we must boost from the emission func- 
tion's frame (i.e. the lab frame) to the pair CM frame, the finite lifetime will 
translate into a non-Gaussian source function through: 



S(r>) 



dr'o 



d 4 RD(R + r/2,p 1 )D{R-r/2,p 2 ) 



(5) 



With the emission function specified, we combine it with the n + ir + kernel to get a 
correlation using Scott Pratt's CRAB code. This kernel includes both the correct 
pair symmetrization and the Coulomb force between the pair. On the left panel of 
Fig.^ we show 2.5 MeV wide slices of the correlation along the sidewards, outwards 
and longitudinal axis. 
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Fig. 1. Left panel: 2.5 MeV slices of the full three-dimensional test correlation. The 
different curves are the test correlation (black), the correlation expanded in Spher- 
ical Harmonics (red), and the imaged, then unimaged, correlation (blue). Right 
panel: terms in the Spherical Harmonic expansion of the test correlation. 

To perform the imaging analysis, we convert the three-dimensional Koonin- 
Pratt Equation in Eq. to a series of one-dimensional equations: 

/>OC 

C em (q) - S eo = An / dr r 2 K e (q,r)Si m (r). (6) 
Jo 

To do this, we simply expand the source function and correlation function in Spher- 
ical Harmonics, expand the kernel in Legendre polynomials, and then apply the 
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Spherical Harmonic Addition Theorem: 
C(q) = ^Yl Yl C z™{q)Y im {q)i S{r) = £ ]T S im (r)Y tm (f), (7) 

£=0 m=-l i=0 m=~l 

and K(q,r) = ^ {21 + l)K t {q, r)P t (q ■ f). (8) 

The right panel in Fig ^ shows the first few terms in the expansion of the test 
problem correlation. 

Having done this, it is natural to ask what these terms mean [[§]. The £ = 
term is the angle averaged correlation and gives us access to i?; n v The £ = 1 
terms give us access to Lednicky offset, in other words which particle was emitted 
first. These terms are absent from like pair correlations due to symmetry. The 
I = 2 terms give us shape information and access to i? , i? s , R\: C20 ~ Rl and 
a±Coo — (C20 ± C22) ~ Rs, Ro (here a± is a constant of order one). Finally, £ — 3 
terms give us access to "boomerang" or triaxial deformation and are absent from 
like pair correlations due to symmetry. Terms with £ > 4 probably will be found in 
future datasets, but are not present in either our test problem or in the real data 
we have analyzed. 

In order to image the test correlation, we must complete converting the one- 
dimensional equations in Eq. JSJ to matrix equations. We do this by representing 
the radial dependence of the source terms in Basis splines: 

N c 

S em (r)=J^S jem B j {r). (9) 
i=i 

Basis splines are piece- wise continuous polynomials and the polynomials that make 
up the splines are patched together at the spline knots. Setting the spline knots is 
both crucial for representing source well and difficult to do. Our method for doing 
this is beyond the scope of this paper and will be discussed in [[7j. 

Inverting the series of matrix equations in Eq. © is an ill-posed problem: the 
kernels are not square and possibly singular, the data is noisy, and error propagation 
is crucial. The practical solution to linear inverse problem, is to find the vector of 
source coefficients that minimize the x 2 '- 

X 2 = {K-S-C) T ■ (A 2 C) _1 -(K-S-C). (10) 

The source that does this is S = A 2 5 • K T ■ (A 2 C)~ 1 • C with the source covariance 
matrix A 2 5 = (K T ■ (A 2 ^" 1 • K)~\ 

3. Data Analysis 

Now we present the results of the imaging analysis on the test problem. In Fig. [21 
we show slices through the imaged source function and the input source function. 
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Clearly we were able to match the height, width and integral of the source function. 
We even were able to begin to resolve the non-Gaussian halo in the longitudinal and 
outwards direction caused by the combination of the finite lifetime of the underlying 
emission function and the lab to pair CM boost. In addition, we can simply read off 
the equivalent Gaussian radii from the source plots and find: i? s = 4.15 ± 0.08 fm, 
R = 5.17 ± 0.12 fm, and R x = 4.53 ± 0.09 fm. We can also read off the height to 
get S(r = 0) = 13 ± 3 x 10 -5 fm~ 3 . If we integrate the source out to the edge of the 
imaged region, we find the fraction of pairs that we can resolve is A = 0.566± 0.025. 
Clearly, we can extract all of the important physics directly from these plots. 
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Fig. 2. Comparison of the imaged source function (blue) and the input source 
function (black). The left set of panels are on a linear scale and the right set of 
panels are on a log scale. 

Finally we apply the imaging to the data from AGS-E859 [Ej. This experiment 
was 14.6 AGeV Si+Au and the data was compiled in the pair CM frame. Fig |3 
shows both the correlation slices and the imaged source slices. Looking at the 
imaged source, we see significant differences between the imaged source and E859's 
best-fit Gaussians. The Gaussian fits and the images mostly agree in the sideward 
and longitudinal directions however the outward direction seems to have a two- 
component structure. If this structure is real, the simplest interpretation is that 
the underlying emission function has an associated lifetime in the longitudinally 
co-moving frame (not the lab frame). Unfortunately, as we see in the correlation 
plots, there may be problems with the data that foil the imaging. There is a large 
band of correlation values substantially below one from q Q = 25 — 75 MeV/c along 
the outwards axis. This is a large unphysical anti-correlation that appears to be an 
experimental artifact [E] to which the Gaussian fits are insensitive. 
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Fig. 3. Left panel: 5 MeV wide slices of the correlations showing the input correla- 
tion (black), input correlation expanded in Spherical Harmonics (red), and imaged 
then unimaged correlation (blue). Right panel: slices through the imaged source 
(blue band) and the best-fit Gaussian from E859 (black). The left set of sub-panels 
are the source on a linear scale and the right set of sub-panels are the source on a 
log scale. 



4. Conclusions 

We have a lot of work to do on this project before it is ready. In addition to 
looking to other datasets, such as the recent PHENIX Run2 data, we are finishing 
the merger of the imaging code with Mike Heffner's CorAL code and Scott Pratt's 
kernel code. In the future we hope to tackle more general test source functions and 
emission functions and to be able to handle un-like particles. 
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